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PREFACE 


The  work  reported  herein  is  sponsored  by  the  Defense  Nuclear  Agency 
under  Contract  No.  DNA  001-77-C-0104  with  Weidlinger  Associates.  It 
constitutes  the  third  phase  of  a study  to  define  the  requirements  for 

applying  dynamic  finite  element  methods  to  analyses  of  MAP  shelter  struc- 

j 

tures.  In  the  current  phase,  the  shelter  analysis  methodology  developed 
in  Phases  1 and  2 of  the  study  is  applied  to  a pretest  simulation  of 
Event  S4  of  the  HAVE  HOST  series.  The  simulation  calculation  is  made  to 
support  the  test  program  and  SAMSO's  closure  design  validation  effort. 

The  authors  would  like  to  thank  Mr.  T.  Kennedy  and  Dr.  C.  Ullrich, 
project  monitors  for  DNA,  for  their  helpful  comments  and  to  SAMSO,  TRW, 
AFWL  and,  in  particular,  the  S4  Test  Working  Group  (S4  TWG)  for  providing 
much  of  the  information  needed  for  the  study. 
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SECTION  I 


INTRODUCTION 

This  is  the  third  phase  of  a study  which  examines  the  application 
of  dynamic  finite  element  methods  to  the  analysis  of  shelter-like  structures 
in  a nuclear  environment.  In  Phases  1 and  2,  the  elastic  deformation  modes 
of  the  shelter  structure,  their  sensitivity  to  variations  in  the  applied  air- 
blast  waveforms  and  the  effect  of  inelasticity  and  computer  modeling 
were  investigated.  An  analysis  methodology  for  performing  dynamic  analy- 
sis of  the  shelter  was  developed.  The  procedure  involves  modeling  the 
shelter-like  structure  and  its  medium  according  to  the  guidelines  estab- 
lished in  the  study,  including  those  characteristics  of  the  air  over- 
pressure which  have  been  found  significant  and  separating  the  analysis 
into  two  distinct  but  complementary  components:  the  short-term  and  long- 
term analyses.  These  and  other  results  from  Phases  1 and  2 have  been 
documented  in  Reference  1. 

The  present  phase  of  the  study  continues  the  work  described  in 
Reference  1,  and  consists  of  applying  the  methodology  developed  to  a parti- 
cular case  study,  namely.  Event  S4  of  the  HAVE  HOST  series.  In  Event  S4 
a half-scale  prototype  shelter  of  the  S4  design  will  be  subjected  to  a HEST 
simulated  attack  environment.  This  simulation  calculation  is  completed 
prior  to  the  test  event,  which  is  scheduled  in  mid-July  1978  at  the  time  this 
report  is  being  prepared. 


The  main  objective  of  the  present  work  is  to  provide  analysis  support 
to  the  test  program  and  to  provide  data  for  the  evaluation  and  validation 
of  the  S4  closure  design.  In  addition,  by  comparing  the  simulation  results 


and  test  data,  when  available,  it  ~ possible  to  obtain  a preliminary 
assessment  of  the  applicability  of  the  analysis  methodology  proposed 
in  Reference  1.  Although  the  methodology  is  general  in  that  it  is  applicable 
to  a class  of  shelter-like  structures  under  various  loading  conditions, 
its  evaluation  based  on  a single  case  study  such  as  the  S4  will  be  re- 
stricted. Of  all  aspects  of  the  methodology,  the  short-term  analysis 
and  headworks  modeling  considerations  can  best  be  measured  in  this  man- 
ner. 

The  remainder  of  this  report  is  organized  as  follows:  A brief  des- 
cription of  the  S4  test  setup  is  given  in  Section  II,  followed  by  a 
description  of  the  simulation  model  in  Section  III.  Because  of  the  vast 
amount  of  simulation  results  obtained,  only  selected  results  pertaining 
to  the  response  phenomenology  are  included  as  Section  IV.  Additional 
simulation  results  which  can  be  used  to  correlate  with  test  data  are  docu- 
mented in  a separate  report  (Reference  11).  Major  findings  and  conclusions 
of  the  study  are  summarized  in  Section  V. 
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SECTION  II 


EVENT  S4 , HAVE  HOST  SERIES 

The  S4  Event  of  the  HAVE  HOST  Test  Program  was  conducted  at 
the  Air  Force  Weapons  Laboratory  (AFWL)  test  site  on  the  Luke  Air  Force 
Range,  Arizona.  A half-size  MAP  shelter  of  the  S4  design  was  tested. 

Unlike  previous  tests  in  the  HAVE  HOST  series  which  utilized  a Dynamic 
Airblast  Simulator  (DABS)  to  provide  a simulation  of  the  airblast  and 
reflected  effects  on  the  door  and  berm  over  the  shelter,  the  S4  Event 
utilizes  the  high  explosive  simulation  technique  (HEST)  to  simulate  the 
local  airblast  acting  at  the  front  face,  headwalls,  and  part  of  the  berm. 
Airblast  induced  ground  motions  were  not  simulated. 

2.1.  S4  Shelter  Configuration 

An  illustration  of  t lie  S4  prototype  is  shown  in  Figure  2-1  (Refer- 
ence 3).  The  shelter  is  buried  with  the  tube  horizontal  and  the  front- 
face  slanting  at  10°  to  the  vertical.  The  details  of  the  headworks  are 
shown  in  Figure  2-2;  it  has  a basically  rectangular  cross  section,  with 
a rectangular  cavity.  Only  the  front-face  and  a very  small  portion  of 
the  top  of  the  headworks  are  exposed.  The  top  is  covered  by  overburden 
and  the  sides  by  two  wingwalls  and  fill.  The  centerline  of  the  rectangular 
headwork  cavity  is  offset  from  the  centerline  of  the  tube. 

The  closure  is  a thick  composite  slab  with  concrete  encased  by  steel 
side  and  back  plates  with  beveled  edges.  One  of  the  main  purposes  of  the 
HAVE  HOST  S4  Event  is  to  provide  test  results  for  the  validation  of  the 
shelter  concept  and,  in  particular,  the  validation  and  evaluation  of  the 
closure/ frame  design. 
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2.2.  HEST  Configuration 

Ttie  HEST  configuration  for  Test  Event  S4  is  shown  in  Figure  2-3 
(Reference  4).  The  two  cavities,  called  the  vertical  and  horizontal 
HEST  cavities,  are  separated  by  approximately  18  inches  of  overburden. 

The  explosive  racks  will  be  preassembled  in  19  modules,  each  2 feet 
wide.  The  modules  will  be  oriented  vertically  on  the  front-face. 

Through  these  two  cavities,  pressure  will  be  applied  to  the  entire 
width  of  the  front-face,  wingwalls,  and  soil  berm  covering  the  headworks. 
Peak  face-on  pressure  is  designed  to  be  4,500  psi,  with  propagation 
velocity  up  the  face  at  37,000  fps.  Peak  top  pressure  is  designed  to 
be  1,400  psi,  with  a velocity  of  7,500  fps  to  simulate  the  actual  airblast 
wavespeed  across  the  top. 

! 

The  HEST  modules  are  constructed  to  yield  pressure  waveforms  to 
match  the  design  waveforms  given  in  Figure  2-4  (Reference  4).  The  peak 
pressure  is  reached  almost  instantaneously  upon  arrival  and  half  of  the 
total  load  impulse  is  delivered  in  the  first  8 msec. 
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Figure  2-3.  HEST  Configuration,  S4  Test 
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Figure  2-4.  WEST  Test  Design  Load  Time  Histories. 


SECTION  III 


SIMULATION  MODEL 

The  finite  element  model  of  the  structure  and  medium  used  in  the 
simulation  analysis  of  the  S4  test  is  basically  the  same  as  that  used 
in  the  main  shelter  study  (Reference  1);  minor  modifications  have  been 
incorporated  to  reflect  differences  between  the  prototype  and  the  test 
configuration  (absence  of  the  floor  slab  and  fill  in  the  tube,  for  example). 
For  completeness'  sake,  a brief  description  of  this  model  will  be  repeated 
herein.  An  approximate  model  for  the  HEST  conf iguration  will  also  be 
discussed  in  this  section. 

One  of  the  key  material  parameters  identified  in  the  main  study  to 
have  significant  effects  on  the  response  of  the  closure  is  the  shear 
strength  of  the  closure  concrete.  Consequently,  in  preparing  the  S4 
simulation  model,  efforts  have  been  spent  in  reviewing  available  data 
on  the  closure  concrete  properties;  the  resultant  model  is  described 
in  Section  3.3. 

3.1.  S4  Shelter  Model 

The  model  assumes  symmetry  about  the  center  vertical  plane  of  the 
tube  and  includes  only  the  left  half  or  hinge  side  of  the  shelter  berm  con- 
figuration. The  finite  element  model  of  the  structure  and  medium  is  shown 
in  Figure  3-1.  It  consists  of  about  18,000  hexahedrons  and,  as  can  be 
seen  from  Figure  3-2  which  shows  the  details  of  the  headworks  model,  special 
attention  has  been  given  to  the: 

(a)  closure  - the  back  and  side  plates  are  modeled  using  two 
elements  across  the  plate-thickness.  The  concrete  portion 
is  modeled  by  three  elements  across  the  thickness. 
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(b)  frame  - the  steel  lining  on  the  inner  frame,  as  well  as 
the  gap  between  the  frame  and  closure  are  modeled.  To 
alleviate  somewhat  the  constraint  on  integration  time 
step  imposed  by  the  miniscule  gap  width,  the  model  incor- 
porates a gap  with  an  artificial  width  of  5 inches  (2.5  inches 
half  scale).  This  change  is  not  expected  to  have  any  effect 
on  the  shelter  response. 

(c)  bearing  ring  - the  steel  bearing  ring  is  modeled  in  the 
refined  fashion  of  the  back  plate.  The  closure,  however, 
is  assumed  rigidly  bonded  to  the  bearing  ring  where  contact 
is  made. 

(d)  transition  section  - the  transition  section  and  the  fore 
portion  of  the  tube  are  sufficiently  refined  so  as  accurately 
to  reproduce  stress  gradients  in  this  region  and  the  trans- 
mission of  the  bending  moment  from  the  headworks  to  the  tube. 

Two  elements  are  used  across  the  thickness  of  the  cylindrical 
wall. 

3.2.  HEST  Load  Model 

The  HEST  configuration  (the  vertical  and  horizontal  cavities,  the  HEST 
overburden,  etc.)  is  replaced  by  the  main  effect  it  is  designed  to  have  on 
the  shelter/berm  configuration,  i.e.,  the  design  loads  as  described  in 
Section  2.2.  Inertial  effects  of  the  HEST  overburden  are  not  simulated  and 
wave  propagation  effects  within  the  HEST  overburden  material  are  approximated 
by  the  application  of  artificial  load-time  histories  on  the  shelter/berm 

surface  which  is  covered  by  the  HEST  overburden  in  the  test  configuration. 


Peak  pressure  is  assumed  to  vary  linearly  in  the  region  between  the  two 


cavities  and  in  the  regions  between  the  cavities  and  the  overburden  limits. 
This  approximation  is  illustrated  in  Figure  3-3  which  shows  the  peak  pres- 
sure distribution  of  the  assumed  loads  on  the  shelter/berm  surface.  The 
times  shown  are  arrival  times  after  first  detonation;  those  corresponding 
to  the  two  HEST  cavities  are  based  on  the  design  detonation  speed  (37,500 
and  7,000  fps,  respectively)  and  those  corresponding  to  HEST  overburden 
^ and  shelter/berm  interface  are  based  on  a wavespeed  of  800  fps  in  the  HEST 

| overburden  material. 

3.3.  Material  Properties 

The  in  situ  soil  material  model  is  identical  to  that  designated  as 
DRY  SAND-1  in  Weidlinger  Associates'  ground  motion  study  (Reference  5), 
with  a density  of  110  pcf  and  an  uniaxial  loading  modulus  corresponding 
to  a wavespeed  of  1,500  fps  and  an  unloading  dililational  wavespeed 
of  3,600  fps.  Its  uniaxial  behavior  is  illustrated  in  Figure  3-4. 

The  backfill  model  is  based  on  data  obtained  from  AFWL  (Reference  6) 
and  corresponds  to  the  so-called  medium  backfill  of  the  data  received. 

f 

The  uniaxial  behavior  of  the  model  is  compared  in  Figure  3-5  with  the  AFWL 
data,  showing  a loading  wavespeed  of  800  fps  and  an  unloading  wavespeed 
of  3,000  fps.  The  dry  weight  density  is  120  pcf. 

The  concrete  is  modeled  as  an  elastic  perfectly  plastic  material,  with 
an  exponential  yield  surface  depending  on  the  mean  compression.  Except 
I for  the  closure,  the  concrete  is  characterized  by  unconfined  compressive 

strength  of  6,000  psi.  The  concrete  in  the  closure,  however,  is  assumed 


to  have  an  unconfined  compressive  strengtli  of  10,000  psi  at  180  days,  in- 


eluding  dynamic  effects  (References  7,  8).  The  yield  surfaces  for  concrete 
in  the  closure  and  in  other  parts  of  the  shelter  are  compared  in  Figure  3-6. 


Note  that  at  a mean  stress  of  10,000  psi  (J^  = 30,000  psi),  the  shear 
strength  of  the  closure  concrete  is  almost  twice  that  of  the  shelter  con- 
crete. 

In  the  main  shelter  study  (Reference  1),  the  shear  strength  of  the 
closure  concrete  is  found  to  have  significant  effect  on  the  response  of 
the  closure  and  this  effect  has  been  discussed  thoroughly.  To  summarize, 
the  closure  behaves  like  a thick  plate  in  flexure  in  the  elastic  range 
(concrete  material  with  infinite  shear  strength)  and  like  a slab  in  punch- 
ing deformation  in  the  inelastic  range  (closure  concrete  witii  limited  shear 
strength,  based  on  6,000  psi  unconfined  strength).  Increasing  the  shear 
strength  by  a factor  of  two  in  the  working  range  of  mean  pressures  will 
tend  to  push  the  closure  behavior  towards  the  elastic  extreme. 

The  structural  steel  in  the  closure  side  and  backplates,  and  in  the 
frame  and  bearing  plate  is  modeled  as  an  elastic  perfectly  plastic  materi- 
al with  a von  Mises  yield  surface  corresponding  to  A36  steel.  No  allowance 
is  made  for  extra  strength  above  the  minimum  due  to  statistical  variation 
or  dynamic  enhancement. 

The  model  parameters  for  the  in  situ  soil,  backfill,  concrete  and  steel 
are  listed  in  Table  3-1.  The  format  is  consistent  with  the  "CAP"  model  des- 
cription, although  only  the  in  situ  soil  and  backfill  are  modeled  as  true 
CAP  materials.  In  t lie  concrete  model  for  parts  of  the  shelter  other  than 
the  closure,  the  tension  cutoff  point  is  varied  based* approximately  on  the 
amount  of  reinforcing  steel  that  exists  in  different  parts  of  the  structure 
(i.e.  , headworks,  frame,  and  tube). 
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3.  A.  Analysis  Program 


The  simulation  calculation  is  performed  using  the  TRANAL  code,  a 
dynamic  inelastic  3D  finite  element  program  developed  by  Weidlinger  Associ- 
ates for  DNA  (Reference  9). 

Since  TRANAL  uses  an  explicit  algorithm  to  integrate  the  equations 
in  time,  the  integration  time  step  used  is  constrained  by  the  mini- 
mum transit  time  across  an  element  for  stability  considerations.  The  sub- 
cycling capability  of  the  code  enables  different  time  steps  to  be  used 
for  different  groups  of  elements  called  zones,  as  long  as  the  stability 
criterion  is  observed  in  each  zone.  In  the  simulation  calculation,  the 
time  step  for  the  zone  containing  the  closure  details  is  necessarily  the 
smallest  due  to  the  small  size  (2.5  inch  full  scale)  of  the  elements  used 
to  model  the  steel  side  and  backplates  and  the  high  wavespeed  of  steel. 

The  computation  time  is  about  50  CDC7600  CPU  seconds  for  0.1  msec  of 
real  time  (full  scale). 

The  simulation  calculation  is  performed  using  the  full  scale  model 
for  the  shelter/berm  configuration.  This  is  done  mainly  for  convenience, 
since  a full  scale  model  of  S4  was  employed  in  the  shelter  study.  The 
HEST  design  loads  are  scaled  appropriately  to  be  compatible  with  the  full 
scale  physical  structure.  The  time  axis  or  impulse  is  multiplied  by  2 
to  go  from  half  scale  to  full  scale.  All  s i mu 1 a t i on  results,  however, 
are  reported  in  the  half-size  or  test  scade  so  they  can  be  compared  readily 
with  test  data. 

The  computation  is  performed  in  installments  of  70  cycles  or  3.5  msec 
(test  scale  time)  each.  Four  installments  have  been  completed  so  far,  and 
simulation  results  with  up  to  14  msec  of  response  time  have  been  obtained. 
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Some  selected  results  will  be  presented  in  the  next  section.  Additional 


results,  mainly  in  the  form  of  stress/strain  invariant  and  time  history 
plots,  are  included  in  Appendix  A;  they  constitute  the  bulk  of  the  data 
provided  for  use  in  the  closure  validation  effort. 
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Figure  '3-6.  Comparison  of  Concrete  Failure  Criteri;: 


MATERIAL  PARAMETERS  FOR  S4  MODEL 


*Vsed  in  conjunction  with  different  tension  cutoffs  (7.2  ksi  to  0.5  ksi)  to  model  various  reenforced 
concrete  in  the  structure,  headwall,  slab  floor  and  grout. 


SECTION  IV 


SIMULATION  RESULTS 


The  response  of  the  S4  test  structure  is,  in  general,  consistent 
with  the  results  reported  in  Reference  1 (Calculation  2A)  for  MAP 
shelters  under  front-on  nuclear  airblast  loading.  The  dominant  deforma- 
tion mode  of  the  structure  is  longitudinal  compression  as  a result  of  the 
loading  on  the  front  face  (the  vertical  HEST).  In  addition,  the  top 
loading  (horizontal  HEST)  pushes  the  headworks  downward.  Detailed  res- 
ponse of  the  test  structure,  however,  differs  from  the  results  of  Refer- 
ence 1.  Some  differences,  such  as  in  the  closure  response,  are  due  to 
differences  in  (closure  concrete)  material  properties  assumed.  Other 
differences,  such  as  in  the  peak  response  parameters,  are  due  to  the  dif- 
ferent loadings  involved.  These  and  other  major  findings  are  discussed 
briefly  in  the  following  paragraphs. 


4.1.  Headworks 

Ttie  deformed  shapes  of  the  headworks  cross-section  in  the  assumed 
plane  of  symmetry  at  several  response  times  are  shown  in  Figure  4-1. 

The  disp lacements  are  exaggerated  for  illustration  purpose.  The  head- 
works  is  pushed  backward  and  downward  by  the  action  of  the  vertical  (front) 
and  horizontal  (top)  HEST  loads.  At  early  times  (before  4 msec  for 
example),  the  top  load  is  not  felt  by  the  structure  and  the  deformation 
mode  is  mainly  that  of  axial  compression.  Slight  downward  motion  is  due 
to  the  10c  incline  of  the  front  face.  The  effect  of  the  top  load  is 
established  8 msec  into  the  simulation  and  the  downward  motion  of  the 
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headworks  can  be  clearlv  identified. 


29 


The  deformed  shapes  of  the  headworks  at  later  times  (12  msec  and 
after)  strongly  suggest  that  the  headworks  moves  as  a rigid-body.  This 
is  not  surprising  since  (as  will  be  shown  later  in  this  section  and  in 


Section  4.2)  much  of  the  closure  and  headworks  remains  elastic.  The 
displacement  profile  at  12  msec  is  typical  of  the  deformed  shape  of  an 
elastic  shelter,  as  described  in  Reference  1. 

The  load  acting  on  the  front  face  induces  high  longitudinal  stresses 
in  the  headworks,  as  shown  in  Figure  4-2.  In  spite  of  these  high  stresses, 
the  headworks  concrete  remains  essentially  elastic  due  to  the  effects  of 
confinement  exerted  by  the  medium  on  the  headworks.  Figure  4-3  is  a 
typical  stress  invariant  plot  for  points  in  the  headworks  and  it  shows 
that  the  headworks  response  is  mostly  in  the  elastic  regime.  Two  excep- 
tions are  the  region  behind  the  bearing  ring  and  the  junction  of  the 
transition  section  and  tube.  The  high  load  on  the  bearing  area  is  a 
reaction  to  the  front  load  acting  on  the  closure,  whereas  the  higli  axial 
stresses  in  the  transition/tube  junction  is  caused  by  an  abrupt  change 
in  load  bearing,  cross-sectional  area  of  the  structure.  They  will  be 
discussed  in  Sections  4.4  and  4.5,  respectively. 

The  longitudinal  stress  in  the  headworks  has  a main  pulse  duration 
of  about  6 msec,  which  is  comparable  to  that  of  the  HEST  load.  At 
points  near  the  transition/tube  junction  (point  D in  Figure  4-2),  the 
pulse  width  is  much  longer  and  the  reason  for  this  will  be  discussed  in 
Section  4.5. 

Since  the  HEST  loads  are  concentrated  in  the  vicinity  of  the  head- 
works,  the  downward  motion  of  t lie  structure  is  most  prominent  there. 
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The  tube  portion  of  the  structure  experiences  very  little  downward 
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motion  and  only  as  a reaction  to  the  downward  motion  of  the  headworks. 

This  is  in  contrast  to  the  response  of  the  shelter  structure  in  a nuclear 
attack  environment.  The  airblast  loading  acting  on  the  berm  induces  down- 
ward ground  (and  hence  structure)  motions  along  the  length  of  the  tube. 

There  is,  of  course,  a time  delay  involved  in  this  action  due  to  the  neces- 
sary transmission  time  through  the  berm  material.  Hence,  at  early  times, 
the  response  of  the  structure  in  a nuclear  environment  can  be  reproduced 
by  that  of  the  test  structure.  The  simulated  and  nuclear  responses  at 
later  times  of  the  structure  as  a whole  and  the  tube  section  in  particular 
may  be  different  due  to  the  absence  of  ground  shock  load  on  top  of  the  tube 
in  the  HEST  test  configuration;  the  extent  of  the  difference  depends  on 
the  backfill  stiffness  (wavespeed)  and,  of  course,  the  intensity  of  the 
load  involved.  For  the  S4  test,  the  ground  shock  will  reach  the  crown  of 
the  fore  portion  of  the  tube  at  about  8 msec  and  the  invert  at  18  msec. 

The  HEST  simulation  response  is,  therefore,  valid  for  at  least  the  first 
8 msec,  or  16  msec  full-scale  time. 

4.2.  Closure 

The  closure  deforms  as  a thick  elastic  slab,  with  maximum  midspan 
deflection  of  the  order  of  an  inch  when  measured  relative  to  the  support. 

This  is  evident  from  Figure  4--4.  The  support  motion,  in  turn,  follows 
the  rigid  body  pattern  of  the  headworks  motion  as  indicated  in  Figure  4-1 
of  Section  4. 1. 

Although  the  overall  deformation  mode  of  the  closure  is  that  of  an 
elastic  slab,  portions  of  the  closure,  notably  the  backplate  and  the  concrete 
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adjacent  to  it,  undergo  inelastic  deformation.  This  can  be  illustrated 
in  Figures  4-5  and  4-6.  The  development  of  inelastic  zones  in  the  clo- 
sure (steel)  backplate  is  illustrated  in  Figure  4-5.  Note  that  the  back- 
plate  material  yields  first  at  midspan,  due  to  the  high  tensile  stresses 
developed  in  bending,  followed  by  yielding  at  the  edges.  Note  also  that 
the  inelastic  activity  starts  at  1.8  msec  or  0.8  msec  after  the  REST 
load  has  attained  its  peak  value  (an  artificial  rise  time  of  1 msec  is 
assumed  in  the  analysis)  and  is  over  2 msec  thereafter. 

The  development  of  inelastic  zones  in  the  closure  concrete  is  il- 
lustrated in  Figure  4-6.  Since  the  concrete  portion  of  the  closure  is 
modeled  by  three  elements  across  its  thickness,  three  profiles  are 
shown  for  each  response  time  of  interest  corresponding,  respectively, 
to  the  front,  mid  and  rear  layers  of  elements  in  the  model.  The  rear 
layer  is  adjacent  to  the  backplate.  Note  that  inelastic  deformation 
in  the  closure  concrete  is  initiated  at  about  1.5  msec  or  0.5  msec 
after  the  load  has  attained  its  peak  magnitude,  and  that  the  initiation 
of  inelastic  deformation  in  the  concrete  precedes  that  in  the  backplate. 
Concrete  at  midspan  fails  in  tension  as  a result  of  flexure;  concrete 
near  the  support  fails  in  shear,  due  to  the  high  shear  stresses  in  the 
slab  and  the  equally  high  reaction  forces  at  the  bearing  plate. 

Note  that  inelastic  deformation  occurs  only  in  concrete  adjacent  to 
the  backplate  (the  rear  layer  in  the  model).  This  finding  is  confirmed 
in  Figure  4-7,  which  shows  the  concrete  failure  criterion  in  stress 
invariant  spaces  (/.!-,',  J^).  Each  data  point  plotted  in  the  figure  cor- 
responds to  the  maximum  response  of  an  element  in  the  closure  concrete 


model.  From  this  figure,  it  is  apparent  that  inelastic  deformation  is 
limited  to  the  rear  concrete  elements.  When  closure  concrete  with  a 
lower  shear  strength  is  used,  as  represented  by  the  dotted  line  in 
Figure  4-7,  it  is  expected  that  inelastic  deformation  will  spread  to 
other  regions  of  the  closure  concrete.  This  is  indeed  the  case  and,  as 
reported  in  Calculation  2A,  Reference  1,  the  closure  deformation  mode 
in  that  case  is  drastically  different  from  that  described  herein. 

The  absolute  velocity/time  histories  of  the  closure  eg  are  given 
in  Figure  4-8.  The  initial  peaks  of  about  500  ips  in  the  longitudinal 
direction  and  180  ips  in  the  vertical  direction  are  mainly  the  response 
of  the  closure  to  the  HKST  load,  because  at  this  early  time  the  headworks 
is  effectively  stationary.  The  portion  of  the  time  histories  after  3 msec 
contains  contribution  from  the  motion  of  the  headworks  (closure  support). 

The  closure  eg  velocity  measured  relative  to  the  support  motion  and 
resolved  into  a component  normal  to  the  backplate  is  given  in  Figure  4-9. 

A Fourier  decomposition  of  this  motion  (Figure  4-10)  shows  a dominant 
period  at  140  Hz.  The  support  or  headworks  motion, on  the  other  hand, 
has  a dominant  period  of  80  Hz  in  both  the  longitudinal  and  vertical 
directions,  as  Figure  4-11  shows.  Dominant  components  at  yet  lower  fre- 
quencies may  exist,  but  because  of  the  short  duration  of  the  calculated 
response  (14  msec)  they  cannot  be  accurately  identified. 

4.3.  Frame 

Minor  inelastic  deformation  occurs  in  the  frame,  notably  in  the  top, 
bottom  and  side  center  gussets,  and  at  the  inside  frame  corners.  Yield- 
ing occurs  at  the  corners  of  the  steel  "ring"  lining  the  inner  surface 
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of  the  frame  (adjacent  to  the  gap  opening)  due  to  the  action  of  the  gap 
pressure  which  tends  to  push  the  frame  outward  and  "round"  the  frame 
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corners.  Two  elements  are  used  to  model  the  lining  across  its  thickness. 

The  inner  steel  element  (adjacent  to  gap)  sustains  high  tension,  whereas 
the  outer  element  (adjacent  to  frame  concrete)  sustains  high  compression 
(see  Figure  4-12).  Concrete  in  the  vicinity  of  the  corners  is  similarly 
affected,  although  to  a lesser  extent.  Inelastic  deformation  in  the 
frame,  however,  is  local  in  that  it  has  no  significant  effect  on  other 
parts  of  the  shelter  so  long  as  the  integrity  of  the  frame  is  not  com- 
promised. 

4.4.  Bearing  Ring 

The  bearing  ring  supports  the  closure  and  is  also  an  integral  part 
of  the  steel  plate  which  forms  the  base  of  the  frame.  The  bearing  load 
time  histories  at  selected  points  in  the  bearing  ring  are  given  in  Fig- 
ure 4-13.  The  bearing  load  is  higher  at  the  bottom  edge  than  at  the  top 
or  side,  a finding  which  is  consistent  with  the  previous  analysis.  Refer- 
ence 1.  Note  that  the  effective  load  duration  is  about  3 or  4 msec,  during 
which  time  the  closure  is  pressed  against  the  bearing  ring  due  to  the  action 
of  the  HEST  load.  The  headworks  is  relatively  stationary  during  this 
time.  Thereafter,  the  closure  and  headworks  move  practically  in  unison 
and  the  bearing  load  is  diminished. 

The  average  peak  bearing  stress  of  25,000  psi  is  about  half  of  the 
static  value  of  54,000  psi  calculated  based  on  the  load/support  area  ratio. 
Furthermore,  the  peak  bearing  load  decreases  rapidly  with  increasing  dis- 
tance away  from  the  closure/bearing  ring  interface.  The  bearing  load  in 
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the  concrete  immediately  behind  the  bearing  ring,  for  instance,  has  a 
peak  of  only  14,000  psi. 

Inelastic  deformation  occurs  in  the  steel  bearing  ring  (Figure  4-14) 
and  the  concrete  immediately  behind  it  (Figure  4-15);  the  extent  of  the 
yielding,  however,  is  slight  and  will  not  impair  the  function  of  the 
bearing  ring. 

4.5.  Transition  Section  and  Tube 

The  high  longitudinal  stresses  induced  in  the  headwords  by  the  front 
load  become  even  higher  in  the  transition  and  tube  regions  because  the 
load  bearing  cross-sectional  area  of  the  structure  decreases.  The  ratio 
of  the  headworks  cross-sectional  area  to  that  of  the  tube,  for  instance, 
is  approximately  5.  Hence,  if  the  strength  of  the  tube  concrete  is  suffici- 
ently high,  a peak  longitudinal  stress  of  about  25,000  psi  can  be  induced 


in  the  tube.  With  a limited  shear  strength,  the  maximum  axial  load  which 
the  tube  can  sustain  depends  on  the  confining  pressure  exerted  by  the  soil 
on  the  tube. 

From  Figure  4-16,  the  average  longitudinal  compression  in  the  tube 
section  a half-tube  diameter  behind  the  headworks  has  a peak  of  8,000  psi 
at  the  crown,  and  10,000  psi  at  the  invert.  Alternately,  the  same  res- 
, ponse  can  be  viewed  in  the  stress  invariant  spaces  (Figure  4-17)  and  it 

is  noted  that  the  peak  longitudinal  compression  is  indeed  limited  by  the 
shear  strength  of  the  concrete  model. 

The  fact  that  the  longitudinal  stress  is  higher  at  the  invert  than 
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at  the  crown  is  significant  in  that  it  deviates  from  the  general  trend 
established  in  Reference  I (Calculations  IB  through  ID,  Calculations  2A 
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through  2C),  although  the  maximum  longitudinal  stress  differential  is 
only  about  2,000  psi.  With  the  floor  slab  and  fill  removed  from  the 
tube,  the  change  in  load  bearing  cross-sectional  area  in  the  test  struc- 
ture is  more  abrupt  at  the  invert  of  the  transition/tube  junction,  (see 
Figure  2-2  of  Section  II)  than  at  the  crown.  This  may  account  for  the 
higher  longitudinal  stress  at  the  invert.  The  HEST  configuration  with 
the  applied  load  restricted  to  the  front  face  and  top  of  the  headworks 
may  also  contribute  to  the  response  described  herein.  This  effect  will 
be  discussed  in  a subsequent  paragraph  of  this  section. 

The  duration  of  significant  longitudinal  stress  in  the  tube  appears 
to  be  longer  than  14  msec,  the  duration  of  calculated  response.  The 
major  cause  of  the  long  duration  is  due  to  yielding  in  the  tube.  The 
longitudinal  load  bearing  strength  of  the  tube  is  limited  because  of 
the  limited  shear  strength  of  the  concrete  and  the  limited  resistance 
provided  by  the  tube  is  not  large  enough  to  arrest  the  momentum  of  the 
headworks  in  a short  time.  It  can  only  slow  down  the  motion  of  the 
headworks  and  has  only  moderate  success  in  doing  so.  Referring  to 
Figure  4-18,  the  mean  headworks  velocity  in  the  longitudinal  direction 
reaches  its  maximum  at  5 msec  and  decelerates  moderately  thereafter. 
Meanwhile,  the  momentum  of  the  headworks  continues  to  carry  it  into  the 
tube,  which  remains  compressed  even  though  the  applied  front  load  is  ef- 
fectively over  in  7 msec. 

From  the  stress  invariant  plot  of  Figure  4-17,  it  is  also  determined 
that  at  about  the  termination  of  the  simulation  calculation,  the  state 
of  stress  in  the  fore  tube  concrete  is  back  into  the  elastic  regime. 
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signaling  the  end  of  inelastic  deformation  (but  not  longitudinal  com- 
pression) at  that  location. 


This  interaction  between  the  headworks  and  tube  can  be  illustrated 
by  the  response  of  a lumped-parameter  model  of  the  structure,  in  which 
the  headworks  is  approximated  by  a mass  connected  to  an  elasto-plastic 
spring  which  represents  the  tube.  The  external  load  is  then  applied  to 
the  mass  directly  and  the  system  parameters  can  be  adjusted  so  as  to 
reproduce  the  overall  headworks/tube  interaction  described  herein. 

The  longitudinal  compression  in  a tube  section  at  about  two-tube 
diameters  behind  the  headworks  is  shown  in  Figure  4-19.  Since  the  com- 
pression is  practically  uniform  throughout  the  cross-section,  only  the 
results  for  points  at  the  crown  are  shown.  The  magnitude  is  about  6,000 
psi,  which  is  attained  upon  wavefront  arrival  and  is  maintained  thereafter. 
Figure  4-20  is  typical  of  the  stress  invariant  plots  obtained  for  points 
in  that  section. 

Comparing  the  stress  invariant  plots  for  the  two  tube  cross-sections 
(Figures  4-17  and  4-20)  brings  out  an  interesting  observation.  Although 
concretes  at  both  locations  deform  ine lastical ly  upon  initial  loading, 
their  responses  differ  thereafter.  For  the  tube  section  near  the  headworks 
(Figure  4-17),  there  is  significant  confinement  pressure  acting  on  the 
lateral  surface  of  the  tube  so  that  higher  mean  pressures  (and  hence 
higher  shear  strengths)  and  higher  axial  loads  can  be  sustained  by  the 
concrete.  The  lateral  confinement  may  come  from  the  engulfment  by  airblast- 
induced  ground  shock,  the  front  load  transmitted  to  the  medium  and  to  a 


lesser  extent  the  restraint  imposed  by  the  massive  headworks  on  the  tube. 


None  of  these  effects  is  apparently  present  at  the  tube  section,  two- 
tube  diameters  behind  the  headworks,  and  consequently,  the  maximum  average 
axial  stress  that  can  be  sustained  is  only  6,000  psi,  which  is  the  unconfined 
strength  of  the  tube  concrete.  It  is  not  known  to  what  extent  this  (lack 
of  confinement)  condition  is  present  in  a nuclear  environment  where  air 
overpressure  loads  act  along  the  length  of  the  berm  and  tube.  The  resul- 
tant soil  confinement  arrising  from  the  top  load  (about  400  psi  for  the 
baseline  nuclear  threat)  depends  on  the  properties  of  the  medium  and  the 
height  of  the  berm  overburden. 

The  strain/time  histories  for  a point  at  the  springline  the  tube 
section  half-tube  diameter  behind  the  headworks  are  shown  in  Figure  4-21, 
and  are  typical  of  the  results  obtained  for  concrete  in  the  fore  portion 
of  the  tube.  The  large  longitudinal  strain  (of  the  order  of  1.52) 
is  due  to  the  longitudinal  compression.  The  large  circumferential  and 
radial  strains  are  a result  of  inelastic  deformation  and  the  dilatancy 
of  t he  concrete  model  assumed;  the  flow  rule  associated  with  the  concrete 
yield  surface  (see  Figure  3-6  of  Section  3.3)  dictates  volumetric  expan- 
sion during  inelastic  deformation. 

Shear  strains  are  in  general  smaller  than  the  normal  strains  and 
are  larger  at  the  springline  (0.5%)  than  at  either  the  crown  or  in- 
vert, where  they  are  almost  negligible. 

Note  that  /ill  strain/time  histories  are  close  to  their  maxima  at 
the  end  of  response  time,  and  because  of  the  predominant  inelastic  defor- 
mation involved,  their  values  at  14  msec  can  be  construed  as  the  accumu- 
lated inelastic  strains. 
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Strains  for  a tube  section  farther  away  from  the  headworks  exhibit 
similar  behavior.  Their  peak  magnitude  decreases  with  increasing  distance 
from  the  headworks.  For  a section  two-tube  diameters  behind  the  headworks, 
for  instance,  the  peak  normal  strain  is  about  O.'IX  in/in  and  the  shear 
strains  are  negligible. 

In  Figure  4-22,  the  displacements  of  a plane  tube  cross-section  are 
illustrated.  The  section  in  question  is  about  one-tube  diameter  behind 
the  headworks,  and  initially,  the  side  or  elevation  view  of  this  section 
appears  as  a straight  vertical  line.  Wien  loaded,  the  higher  longitudinal 
stress  at  the  invert  is  manifested  as  a larger  longitudinal  displacement 
at  that  location  than,  sav,  at  the  crown.  II  straight  lines  are 
drawn  through  the  data  points  to  represent  gross  motion  of  the  tube  sec- 
tion had  it  remained  plane,  these  lines  will  rotate  in  the  clockwise 
direction  signifying  bending  of  the  tube  in  that  direction.  In  the  actual 
deformed  states,  as  shown  in  Figure  4-22,  most  of  the  non-plane  deforma- 
tion occurs  in  the  invert  region. 

Note  that  at  early  times,  the  crown  of  the  section  is  actually  dis- 
placed upward,  as  a result  of  the  longitudinal  compression  and  dilatancv 
of  concrete.  In  fact,  the  tube  section  expands  radially  outward  as  the 
tube  is  compressed  in  the  longitudinal  direction. 

I'he  deformed  profiles  of  a tube  section  about  two-tube  diameters  be- 
hind the  headworks  is  given  in  Figure  4-2  1.  Hie  deformation  is  quite 
uniform  throughout  the  section,  consistent  with  the  uniform  distribution 
ol  longitudinal  stress  observed  in  Figure  4-19.  The  displacement  is  less 


than  that  for  a section  closer  to  the  headworks  (see  Figure  4-22),  also 


consistent  with  the  lower  longitudinal  stress  recorded  at  that  location 
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(see  Figure  4-19).  There  is  also  less  expansion  of  the  section  in  the 
radial  direction. 

Circumferential/time  histories  for  points  in  a tube  section  half- 
tube diameter  behind  the  headworks  are  shown  in  Figure  4-24.  The  mean 
hoop  compression  is  initially  higher  at  the  invert  than  at  the  crown  or 
springline.  This  is  a reaction  to  the  longitudinal  compression  since,  as 
noted  in  Figure  4-16,  higher  longitudinal  stress  is  induced  at  the  invert 
than  at  the  crown  or  springline.  Comparing  the  circumferential  stresses 
at  the  inner  and  outer  fibres  at  the  crown  of  the  tube  section  (Figure  4-24a) 
indicates  that  the  ring  segment  at  that  location  undergoes  initially  nega- 
tive (convex)  bending  and  later  positive  (concave)  bending  in  the  plane 
of  the  section.  Applying  the  same  reasoning  to  the  invert  and  springline 
(Figures  4-24h,c),  the  ring  segments  at  the  invert  and  springline  are  found 
to  sustain  inplane  convex  bending  only. 

The  inplane  deformation  can  also  be  illustrated  by  Figure  4-25,  in 
which  are  plotted  the  deformed  shapes  of  the  tube  cross-section  in  question. 

At  early  times  (before  10  msec),  the  inplane  deformation  is  mainly  radial 
as  the  tube  expands  outward  when  compressed  longitudinally.  The  inplane 
bending  patterns  at  the  crown,  invert  and  springline  described  in  the  pre- 
ceding paragraph  are  also  apparent.  They  consitute  the  n = 4 deformation 
mode  of  a ring,  and  this  deformation  mode  of  relatively  high  order  is 
probably  induced  by  the  structural  constraints  imposed  by  the  bulky  headworks  on 
the  tube.  After  10  msec,  the  downward  motion  of  the  headworks  (and  top 
MKST  load)  is  established  and  its  effect  is  to  pull  the  tube  downward  also. 
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The  inplane  deformat  ion  mode  is  now  dominated  hy  the  oval  ling  or  n = 3 
mode  except  near  the  invert  where  the  stiff  in  situ  soil  in  contact  with 
the  structure  constrains  it  from  following  the  oval  ling  pattern. 

The  inplane  deformation  pattern  for  a tube  section  farther  down  the 
tube  (two-tube  diameters  behind  the  headworks)  is  illustrated  in  Figure  4-2h. 
Because  of  the  large  distance  away  from  the  headworks,  the  effect  of  the 
structural  constraint  of  the  headworks  on  the  early  time  deformation  mode 
is  not  as  prominent  as  that  shown  in  Figure  4-25.  The  late  time  deforma- 
tion pattern  is  similar,  however,  to  that  of  a section  near  the  headworks, 
although  the  magnitude  of  deformation  is  proportionally  less.  (Note  that 
the  displacement  scale  in  Figure  4-25  is  about  twice  that  of  Figure  4-26.) 

4.6.  Soi I /St ructure  Interface  Stresses 

The  distribution  of  soil  pressure  on  t lie  bottom  side  of  the  headworks 
is  illustrated  in  Figure  4-27;  the  peak  value  ranges  from  400  to  600  psi. 

The  term  pressure  is  used  in  this  section  to  denote  the  stress  component 
normal  to  the  soil/structure  interface.  Note  that  the  initial  peak  for 
the  soil  pressure  at  point  A at  3.5  msec  is  due  to  the  REST  induced  ground 
shock,  whereas  the  later  peak  at  10  msec  is  probably  due  to  the  downward 
movement  of  the  headworks.  Soil  resistance  to  structural  motion  increases 
with  distance  from  the  front  (from  A to  D). 

Soil  pressures  acting  on  top  of  the  headworks  are  given  in  Figure  4-28. 
Since  they  originate  from  the  top  REST  load,  the  magnitude  and  the  character- 
istics of  t lie  soil  pressure  at  a particular  point  on  top  of  the  headworks  de- 
pend significantly  on  the  thickness  of  the  berm  above  that  location.  For  a 
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point  tinder  a shallow  layer  of  backfill,  such  as  point  A of  Figure  4-28, 
the  induced  soil  pressure  is  high  and  the  pulse  characteristics  are  similar 
to  the  HEST  load  (spurious  behavior  after  the  initial  peak  is  due  partly  to 
reflections  at  the  soi /st ructure  interface  and  the  berm  face  surface,  and 
partly  to  numerical  integration).  At  point  B of  Figure  4-28,  the  effect  on 
the  pulse  shape  of  traveling  through  five  feet  of  soil  is  conspicuous;  the 
peak  soil  stress  decreases  to  1100  psi  and  the  pulse  shape  is  significantly 
modified.  These  simulation  results  include  some  "numerical  filtering"  due 
to  the  finite  dimension  of  the  elements  used.  However,  they  provide  a good 
measure  of  the  load  (impulse)  responsible  for  the  downward  motion  of  the 
headworks,  since  that  motion  has  a basic  frequency  (80  Hz  from  Figure  4-11) 
well  within  the  capability  of  the  model. 


Soil  pressures  acting  on  the  sides  of  the  headworks  are  illustrated 
in  Figures  4-29a  through  4-29d,  at  various  distances  from  the  point  and  at 
two  different  depths.  Peak  pressure  is  about  2,000  psi,  and  is  due  mainly 
to  stress  waves  induced  in  the  medium  by  the  front  HEST  load.  The  peak 
pressure  decreases  with  distance  from  the  front  face  and  is  only  800  psi 
at  the  transition/tube  junction.  The  pulse  shape,  however,  remains  relative- 
ly unchanged  with  distance. 

Soil  pressures  acting  on  the  tube  at  early  times  (before  14  msec)  are 
passive  pressures,  i.e.,  they  are  exclusively  due  to  the  actions  of  the 
front-loaded  structure  against  the  soil.  Their  distribution  along  the  tube 
is  illustrated  in  Figure  4-30  for  points  at  the  crown,  in  Figure  4-31  for 
points  at  the  invert  and  in  Figure  4-32  for  points  at  the  springline.  The 
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passive  pressure  is  higher  at  the  front  portion  of  the  tube  and  becomes 
negligible  at  about  two-tube  diameters  behind  the  headworks  and  thereafter 
This  lack  of  confinement  at  mid-span  confirmed,  the  explanation  given  in 
Section  4.5  for  the  low  longitudinal  load  sustained  by  the  tube  at  that 


location. 


Deformation  Mode  of  Headworks,  Plane  of  Symmetry. 


Figure  4-2a.  Longitudinal  Stress/Time  Histories  for  Points  Along  Upper  Headworks. 


Figure  4-3.  Typical  Stress  Invariant  Plot  for  Concrete  in  Headworks. 


Inelastic  Zones  in  Closure  Concrete,  Rise  Time  of  Load  is  1 Msec  (Continued). 


Figure  4-7.  Distribution  of  Peak  Closure  Concrete  Stresses  in  Stress  Invariant  Space 


Velocity /Time  Histories,  Closure  Center  of  Gravity 


Relative  to  Support  Motion. 


Figure  4-lla.  Fourier  Decomposition  of  Average  Closure  Support  Motion,  Vertical  Component. 


Figure  4-12a.  Stress  Invariant  Plot  for  Steel  in  Upper  Frame  Corner,  Inner  Lining 


Figure  4-12b.  Stress  Invariant  Plot  for  Steel  in  Upper  Frame  Corner,  Outer  Lining 


Figure  4-14.  Typical  Stress  Invariant  Plot  for  Steel  in  Bearing  Ring. 
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Shear  Strain/Time  Histories  for  a Point  at  Springline  of  Tube  Section  Half-Tube  Diameter 
Behind  Headworks  . 
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Figure  4-22. 


Displaced  Profiles  of  Tube  Section,  One-Tube  Diameter 
Behind  Handworks. 


Displaced  Profiles  of  Tube  Spot  Ion  Two-Tube  Diameters  Rebind  'leadworks 


Figure  4-24.1.  Circumferential  Stress/Time  Histories  for  Points  at  Crown  of  Tube  Section 
Half-Tube  Diameter  Behind  Headworks. 


Half-Tube  Diameter  Behind  Headworks 


ube  Diameter  Behind  Headworks. 


5.  Deformed  Shapes  of  Tube  Cross-Section  One-Tube  Diameter 
behind  tleadworks. 
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SECTION  V 


SUMMARY 

The  dominant  deformation  mode  of  the  test  structure  is  longitudinal 
compression,  with  some  downward  movement  of  the  headworks  at  later  times 
(after  8 msec).  Despite  high  stresses,  the  bulk  of  the  headworks  and  clo- 
sure remains  elastic;  minor  inelastic  deformation  occurs  in  the  closure 
backplate  and  its  adjacent  concrete  in  the  closure,  the  center  gussets  and 
inside  corners  of  the  frame,  the  bearing  ring  and  the  concrete  immediately 
behind  it. 

Significant  inelastic  deformation  occurs  in  the  tube  section;  the  de- 
formation is  more  severe  for  portions  of  the  tube  closer  behind  the  headworks. 
The  tube  section  about  one-half  tube  diameter  behind  the  headworks  sustains 
longitudinal  stresses  of  the  order  of  10,000  psi  and  (plastic)  strains  of 

the  order  of  1.5%.  The  tube  section  about  two  tube  diameters  behind  the 
headworks  sustains  longitudinal  stresses  of  the  order  of  6,000  psi  and 
strains  of  the  order  of  0.3%. 

The  S4  REST  load  configuration  does  not  provide  ground  shock  loading 
on  the  top  berm  along  the  length  of  the  tube.  This  constraint  has  two  impor- 
tant impacts  on  the  test  structure  response. 

. Downward  motion  of  the  structure  is  most  prominent  at  the  headworks; 
the  tube  portion  experiences  very  little  downward  motion  and  only  as 
a reaction  to  that  of  the  headworks.  The  downward  headworks  motion 
dominates  the  late  time  (after  10  msec),  in-plane  deformation  pattern 
of  the  tube  cross-sections,  which  is  mainly  the  ovalling  mode  except 
in  the  vicinity  of  the  invert  where  the  stiff  in  situ  soil  in  contact 
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with  the  tube  cross-section  constrains  it  from  following  the 
ovalling  pattern.  The  average  downward  tube  displacement  de- 
creases rapidly  with  distance  from  the  headworks.  At  1A  msec 
it  is  about  an  inch  for  a tube  section  one  tube  diameter  behind 
the  headworks,  and  about  0.3  inches  for  a section  another  tube 
diameter  farther  back. 

. Portions  of  the  tube  farther  than  one  diameter  behind  the  head- 
works  are  effectively  unconfined,  i.e.,  the  active  earth  pressure 
is  null,  and  the  tube  concrete  yields  at  an  average  longitudinal 
compression  of  6,000  psi,  its  assumed  unconfined  strength. 

The  closure  concrete  model  used  in  the  simulation  has  a yield  strength 
much  higher  than  that  used  in  Reference  1 (see  Figure  3-6).  This  difference 
has  two  major  effects. 

. The  closure  in  EVENT  SA  responds  basically  as  an  elastic  slab, 
with  nominal  mid-span  deflections  of  an  inch.  This  deformation 
mode  is  in  contrast  to  the  inelastic  punching  deformation  observed 
in  Reference  1 for  a closure  with  low  strength  concrete.  This  change 
in  the  closure  response  is,  in  part,  influenced  by  the  lower  HEST 
load  simulated  in  the  test  (A, 500  psi  for  SA  test  versus  5,300  psi 
used  in  Reference  1),  but  the  effect  of  the  concrete  properties 
appears  to  be  more  significant  than  that  of  the  peak  loading. 

. By  remaining  elastic,  the  SA  closure  is  able  to  transmit  most  of  the 
front  load  impulse  to  the  headworks,  thereby  giving  the  latter  a large 
momentum.  This  momentum,  when  combined  with  the  lack  of  active  earth 
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pressure  at  the  tube  portion  of  the  test  structure,  is  responsible 
for  the  significant  inelastic  deformation  observed  in  the  tube. 

The  absence  of  the  floor  slab  and  fill  inside  the  tube  of  the  test 
structure  results  in  an  abrupt  change  in  the  load  bearing  cross-sectional 
area  at  the  junction  of  the  transition  section  and  the  tube,  this  change 
being  more  abrupt  at  the  invert  than  elsewhere  at  the  tube  cross-section. 
Consequently,  the  longitudinal  stress  is  amplified  as  it  propagates  from 
the  transition  section  to  the  tube,  more  so  at  the  invert  than  at  the  crown 
or  springline.  This  creates  a slight  but  distinct  longitudinal  bending  moment 
which  tends  to  torque  the  tube  in  a sense  (convex)  which  is  opposite  to  the 
case  where  the  fill  and  slab  are  present.  This  torque,  however,  diminishes 
rapidly  with  distance  along  the  tube  and  is  negligible  at  a tube  section 
two-tube  diameters  behind  the  headworks. 
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Air  Force  Systems  Command 
ATTN:  XRTB 

Strategic  Air  Command 
Department  of  the  Air  Force 
ATTN:  XPFS 

ATTN:  NRI-STINFO  Library 

Assistant  Chief  of  Staff,  Intelligence 
Department  of  the  Air  Force 
ATTN:  IN 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS 


Civil  Systems  Inc. 

ATTN:  J.  Bratton 


General  Electric  Company-TEMPO 
ATTN:  DAS  I AC 

University  of  Illinois 
Consulting  Services 

ATTN:  N.  Newmark 

Karagozian  and  Case 

ATTN:  J.  Karagozian 

Merritt  CASES,  Inc. 

ATTN:  J.  Merritt 

University  of  New  Mexico 

Electrical  Engineering  & Computer  Science  Dept. 
ATTN:  G.  Lane 

Physics  International  Co. 

ATTN:  F.  Sauer 

R & D Associates 

ATTN:  C.  MacDonald 
ATTN:  J.  Lewis 

Rand  Corp 

ATTN:  C.  Mow 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

Science  Applications,  Inc. 

ATTN:  H.  Wilson 

Science  Appl ications , Inc. 

ATTN:  J.  Dishon 

Science  Applications,  Inc. 

ATTN:  D.  Maxwell 

SRI  International 

ATTN:  G.  Abrahamson 

Systems,  Science  & Software,  Inc. 

ATTN:  D.  Grine 

Systems,  Science  & Software,  Inc. 

ATTN:  C.  Hastings 

Tech  Reps 

ATTN:  R.  Holmes 

TRW  Defense  & Space  Sys.  Group 
2 cy  ATTN:  P.  Dai 

Weidl inger  Assoc . , Consulting  Engineers 
ATTN:  J.  Isenberg 
ATTN:  F.  Wong 

Boeing  Co. 

ATTN:  R.  Holmes 
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